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C^^ \ The E° — A mixing angle in isospin-asymmetric nuclear medium is inves- 

tigated by using QCD sum rules. From the general consideration of the in- 

^P \ medium baryonic correlations, in- medium baryon mixings are shown to have 

several Lorentz structures such as the scalar mixing angle 6 and the vector mix- 
ing angle 6^ . This causes a difference between the particle mixing (= 6^ + 9^) 



D \ and the anti-particle mixing [= 6^ — 0^). From the finite energy sum rules 



for the S° — A mixing, we find that the in-medium part of the mixing angle has 
a relation 9^^^ ~ ~^Mcd i'^ the isospin-asymmetric medium. This implies that 
the medium affects mainly the anti-particle mixing. From the Borel sum rules, 
we obtain \9 — 9q\ ~ 0.39 | {pn — Pp)\/ po with 9^, pn, Pp and po being the vacuum 
mixing angle, the neutron density, the proton density and the normal nuclear 
matter density respectively. 
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1 INTRODUCTION 

The SU{2) isospin symmetry is slightly broken in the hadronic world. Examples of 
this symmetry breaking are the mass splittings within a same isospin multiplet {p — n, 
S^ — S" and vr* — vr"), the particle mixing among different isospin multiplets (vr" — i], 
71° — T]', p° — uj and S" — A), and the nuclear force in the ^5*0 channel {Vpp ^ Vpn 7^ Vnn) 

ii- 

The isospin symmetry breaking has two different sources: (i) the electromag- 
netic (EM) effect due to the electric-charge difference between u and d (e„ 7^ e^), 
and (ii) the quark-mass difference between u and d {rriu 7^ rrid)- The latter effect can 
be evaluated from the mass term TiSas? of ^^e QCD Hamiltonian density for light 
flavors, 

"^SS = 2("^« + ^d)(uu + dd) + -{rrid - mu)(dd-uu) + m^ss, (1.1) 

where (m^ — mu)/{rnu + rrid) ~ 0.29 0]. Ti^ass i^ known to be more important than 
the EM effect for the p — n mass difference, the p° —uj mixing, and the S" — A mixing 

i- 

The QCD sum rule is a useful method to evaluate the magnitude of the 

isospin symmetry breaking with non-perturbative QCD dynamics. It has been applied 

for the isospin mass splittings in octet baryons 1§, 0, |§], the p° — u mixing [^, ^, the 



n° — 7] mixing |10| and the S" — A mixing [|Tl|, |T^. For example, the E° — A mixing 



angle in the vacuum defined by 






EJQCD 



° Mso - Ma 



a: 

1.2) 



is evaluated as |6'o| = 1.4 x 10 ^ |]T^ and 7 x 10 ^ ||TT| in QCD sum rules. This value 



is comparable to the other estimate |^u| :ii 1 x 10^^ in the naive quark model |jT3[ and 
in the chiral perturbation theory ^, Q. 

In this paper, we will consider the isospin-asymmetric nuclear medium where 
the difference between the neutron density (p„) and the proton density (pp) becomes 
an extra source of the isospin symmetry breaking. In particular, we study how this 



new source affects the S° — A mixing. The in-medium QCD sum rule |T^, |T5|, ^ is a 
suitable method for this purpose, since we can treat the isospin-asymmetric medium 
as a background field acting on the S" — A correlation through the operator product 
expansion. Also, it allows us to investigate the response of the mixing angle under 
the variation of the magnitude of isospin-asymmetry. 

We should mention here that the p° — u mixing in isospin-asymmetric medium 
has been recently studied in Jl^. A major difference between the meson-mixing 



treated in WTj and the baryon-mixing in the present paper hes in the fact that the 



latter can have several Lorentz structures (such as scalar and vector mixing angles) 
in the medium because of the spinor structure of the baryon fields. This will cause 
an interesting difference between the particle mixing (S" — A) and the anti-particle 
mixing (S" — A). In the former (latter), the scalar mixing and the vector mixing 
act in destructive (constructive) manner. A close analogy of this phenomenon is the 
scalar and vector self-energies of the nucleon (anti-nucleon) in the symmetric nuclear 



medium, where scalar and vector act in destructive (constructive) way fl^ . 

The organization of this paper is as follows. In SecH, we analyze the general 
structure of in-medium correlation functions of spin ^ baryonic currents with and 
without the mixing. Dispersion relations satisfied by the correlation functions are 
also written down after decomposing them into even and odd parts with respects to 
the frequency uj of the currents. In Sec.|^, we introduce a generalized mixing matrix 
in the spinor space and make physical interpretation of the scalar and vector mixing 
angles. In Sec.§, we carry out the operator product expansion (OPE) of the mixed 
correlation function in the S°— A channel. The Lorentz-tensor and isospin-asymmetric 
operators are kept in OPE since they have non-vanishing expectation values in the 
isospin-asymmetric medium. The in-medium condensates which appear in OPE are 
also evaluated in the low density expansion in this section. Since the in-medium 
expectation values of the isospin-asymmetric operators beyond dimension 4 are hard 
to be determined at present, we limit ourselves to the OPE up to dimension 4. In 
Sec.|^, we construct the finite energy sum rules [|1^] and the Borel sum rules using 
the results in previous sections. Then we extract a qualitative result from the finite 
energy sum rule. In Sec.^, to reduce the uncertainties due to the absence of higher 
dimensional operators in OPE, we examine the reliability of the sum rules constructed 
in Sec.|^ from the point of view of the consistency among different sum rules. Then 
we evaluate the E" — A mixing angles numerically. Sec.|^ is devoted to summary and 
concluding remarks. 

2 GENERAL PROPERTIES OF THE CORRE- 
LATION FUNCTIONS 

In this section we examine the spinor structures of the diagonal and off-diagonal 
correlation functions of spin ^ baryonic currents. We will also derive the dispersion 
relation for each spinor component of the correlation functions. 



2.1 Spinor structure 

Let us start with the following two-point functions: 

n^(p|g) = ^|rf1re^P^ (g|T[r/^ (x) 77^(0)] |g), (2.1) 

U^{p\q) = tJd'xe'^^q\R[vA{x)rjM]\q), (2-2) 

where T and R denote time-ordered and retarded products respectively, and riA{B) (x) 
is an interpolating operator for the baryon A{B). If A is different from B, the 
correlations describe the particle mixings. |g) is a state vector with four-momentum 
qf^. Later, this state will be identified with the isospin- asymmetric nuclear medium 
to investigate the S" — A mixing with A = A and S = S". 

These correlation functions have the following spectral representations, 



1 f°° 
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(2.3) 
(2.4) 



with p^ = {p° ,p) and p'^ = {p'", p). p^^ {p \ q) and p^,^ {p \ q) are the spectral functions 
defined by 



Pa/3 {P I Q) 
Pc^P {P I Q) 



d^e^^^(q{Mx)}AVBm, 
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dxe 



{VBmp{VA{x)]^ 



(2.5) 
(2.6) 



where spinor indices (a,/3) are explicitly written. 

To make the following discussion concise, let us introduce a linear combination 
of the spectral functions with real parameters a and h as 



and define 11^ as 



A(p| q) = ap{p\q) + hp {p\ q) 

.,0 A(p'k) 



1 /-oo / 
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Then the time-ordered and retarded correlations can be written as 



a=0 

6=1/2 



Ii^{p q) = 


-- U_ {p q) ,^1/2 + n+ (p q) 

fe-O 


n^(p q) = 


= n (n g) 

- ^^ ^'' a=b=l/2. 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



Let us first consider the spinor structure of the spectral function (|2.7| ). A has 
a 4 X 4 spinor structure and can be expanded in terms of a complete set of Dirac 



matrices. The Lorentz covariance restricts the general form of A as 

A{p\q) = A^ + AV + A^> + A^'^ + A^':^7^ + A^'^7^ 

+A^H^^m^pV) + A^^(a>V), (2.11) 

where we define the coefficients A' = A\p'^,p-q,q'^) for / = S, P,Vi,V2, Ai, A2,Ti,T2 
and a^^ = \e^yai3(^°'^ = —io'fj.ul^ with a convention e^^'^^ = 1. 

Parity and time-reversal properties further restrict the spinor structure of A. 
Since the baryonic currents 77^(2^) and t]b{x) have the same transformation properties 
with elementary Dirac fields under parity {V) and time-reversal (T) transformations, 
we have 

Vr]A{x)V-^ = Pr]A{S), Vr]A{x)V-^ = r]A{x)p-\ (2.12) 

Tr]A{x)T-' = Tr]A{-S), rr]A{x)r~'=r]A{-x)T-\ (2.13) 

where x'^ = {x°, x) and x'^ = {x'\ —x). P (T) is a 4 x 4 matrix in the spinor space 
for the parity (time-reversal) transformation. We assume the same transformation 
matrices P and T for riB{x). Under the Hermitian conjugate, P and T, the Dirac 
matrices transform as 

Y [{1,7', 7m, 7/.7', «^/.., (^lu}]\° = {1, -7', 7/., 7/.7', -^^M-, ^'.} , (2.14) 

P{l,7^ 7^, 7^7^ ia^„ aj,}p-i = {1, -7^ %, -%7^ z?^,, -aj,} , (2.15) 

T{l,7^ 7^, 7^7^ ^a^„ aJ^JT-^]* = {l, 7^ 7^, 7^7^ ^?^„ SrjJ , (2.16) 

where 7^" = (7^,7), 7^ = (7°, -7), ^f^u = ^[lf,,lu\ and a^^ = -ia^^j^. The state 
vector |g) is assumed to have the property, 

V\q) = \q}, T\q) = \q), (2.17) 

with q^ = {q°, q) and q^ = {q°, —q). 

Because of the transformation properties of the baryonic current and the state 
vector shown above, the spectral function satisfies the following relations 

A(p|g)U^5 = 7"{A(j5k)}S°, 

A(p|g) = PA{p\q)p-\ (2.18) 

A{p\q) = {TA{p\q)T-'y, 

where A^-^B stands for the exchange of A and B. 

Owing to Eq. (|2.18|) and Eqs.( p.l4| - p.l6D , eight independent functions in Eq.( |2.11|) 



reduces to four functions such as 

A(p I g) = A^ + A^'p + A^^q + A^i {la^.p^q''), (2.19) 



where 



A' 



A^B 



A' {l = S,VuV2), A^i 



A*^B 



-A^\ 



(2.20) 



A'* = A' il = S,V^,V2,T,). 
Eqs.( p.l9 ) and (|2.20|) also imply that the correlation function 11 defined in Eq.(p.8| 



has a form 



n^ (p I g) = n J + Ul'p + Ul'q + Ul' {la^^p'^q") 



where Hi- = II' (p^, p-g, q^) for I = S,Vi, V2, Ti and 



n' 



A^B 



n' {i = s,v,,V2), ul^ 



A^B 



-ul\ 



(2.21) 



(2.22) 



Because of Eqs. (p.9|) and ( p.lO|) , the time-ordered and retarded correlation functions 
have the same decomposition as Eq.( |2.21| ). 

Note that, for A = B, our results are fully consistent with the previous analysis 
in Ref.||20|. In particular, the tensor terms A^^ and 11^^ vanish in this case. For 
A ^ B, our results are new. The tensor terms do not vanish in this case unless p'^ 
and q^ satisfy special conditions. 



2.2 Dispersion relations 



Eq. (|2.21|) enables us to decompose the dispersion relation (p.8|) into independent 
structures. In this subsection we will work in the rest frame of the state vector \q) 
{q^ = (g°, 0)), since it is sufficient for later applications. For notational simplicity, we 
will omit the argument q in 11^, whenever we consider the rest frame of |g). 
Under this simplification, Eq.( |2.2lD becomes 



n^ ip) = u'+ n^7° - (n^^ - ul^q"!") {p-i) 



where we have introduced 



n 



V 



Yi^ p° + n^2 q°. 



Then the dispersion relation for I = S,V reads 

Ren^(a;,|p|) 
where 



±-p/"do;'^"^^^ (^'> IpI) 



vr 



UJ' — UJ 



ImUL {uj,p) = ±— lim [RcII' (u + ie, \p\) - ReUL {cu - ie, \p\) 



'2i e-^+o 



(2.23) 
(2.24) 

(2.25) 
(2.26) 



with UJ = p" and P stands for the principal value integral. As far as p 7^ 0, the same 
dispersion relation holds for / = Vi,Ti. 



Next let us decompose the correlation functions to even and odd parts under 
the transformation uo ^^ —uo [p]| : 



n^ {u, \p\) = Yi'i^^ (u\ \p\) + cu n^- ' yu^, ipi 



r'(0) (. ,2 



(2.27) 



Then Eq. ( p.25| ) reduces to a formula which relates the even (odd) part of RerT-j- with 
the odd (even) part of Imfl-j-: 



no) (, 



1 /-oo Imn'p^^^(sMp|)v^ 



s' — s 



(2.28) 



Renli^^ (s, Ipl) = ±-Frds' 

71 Jo 

Renr(^.li>l) = ±ip/V-""""'--l''l>^^. (2.29) 

71 Jo S' — S 



Also, Eq.( p.26D reduces to 



imn 



HO) 



Imn 



1{E) 



2i e^+O 

1 



W(E) 



s, \p\) = ±— lim Ren'^^''^ {s + ie, \p\) - Ren^^""^ {s-ie, \p\) /v^, (2.30) 



1{E) 



T«(0) 



s, \p\) = ±- lim Ren^^""^ {s + ie, \p\) - ReK^''> {s-ie, \p\) ^s, (2.31) 



l{0) 



'2i e^+o 
for I = S,V, V\, Ti where s = uP' . 

2.3 Dispersion relations between Rellrp and Imllj^ 



The retarded correlation function defined in Eq. (|2.1CI| ) satisfies the same dispersion 
relations (|2:25D , ( ^12^ ) and dOOl) . Also, Eqs.(|2j) and ^^ imply that the real part 
of IIj- and that of Vl^ are equal for each spinor component: 

Ren^ = Ren^ (/ = S,Vx, V2, Ti,V) . (2.32) 

Therefore, the dispersion relations in the rest frame of the state vector |g) reads 

Ren^ {uj, \p\) = -P du' ^^^1^^ (2.33) 

71 J -00 Uj' — UJ 

Imn^(^,p) = 



— lim 

2i e^+o 



Rellr {oj + ie, \p\) - Rell^ {u - ie, \p\) 
By decomposing the above to the even and odd parts, one finds 

W(o) 



TT JO S' — S 

ReUr> (., \P\) = ip/-,.,i!£HZVMZ^ 

^ ^ ' '^ 7T Jo S' -S 



(2.34) 

(2.35) 
(2.36) 



ReU^r}^'' {s+ie, p ) - 


-Re4(^)(s-^e,|p)]/v^, 


(2.37) 


Ren^^'^) {s+te, p\) - 


- Renfr^°^ (s-^e, p|) 


v^, 


(2.38) 



and 

Imlli (s, Ipl) = — lim 

Imn^^^^ (s, IpI) = ^£m 

/ = 5, y, \4, Ti with n^^^ = v^T\Rf + n^'ij g". 

In the sum rule analysis in later sections, Relly is evaluated by the operator 
product expansion, while phenomenological ansatzes are made for Imll^. 

3 STRUCTURE OF BARYON MIXING 

In this section, we discuss the general structure of the baryon mixing. 

3.1 Definition of the mixing angle 

Let us consider a mass matrix for particles A and B which have definite quantum 
numbers but have different masses, 

m4^!'^\h{\A) |B))=f, ''\ ^ ^m\^)\^ ,3,1) 

\{B\ ) V I / I / y y{A\E\B) Mb j 

where H is the Hamiltonian of the system and \A) and \B) are one-particle states 
normalized as (A|A) = {B\B) = 1. The diagonal matrix element is equal to the 'mass' 
of the symmetric state, namely M4 = {A \H\ A) and Mb = {B \H\ B). We choose the 
relative phase of the states \A) and \B) to be {B \H\ A) = {B \H\ A)* = {A \H\ B). 
The physical states |^)phys and \B)^,^ are represented as a linear combination of \A) 
and |i?) by using the mixing angle 9, 

l^)phys = |A)cos^+|i?)sin^, (3.2) 

l^)phys = |5)cos^-|A)sin^. (3.3) 

The mass matrix in terms of the physical states reads 

_/ M-AM'/2 {A\H\ B) COS29 + AM sm2e/2\ 

^^"""'y {A\H\B) COS29 + AM sm2e/2 M + AM' /2 )^ ' ' 

where M={Ma + Mb) /2, AM = Mb- Ma iind AM' = AM cos29 -2 {A\H\ B) sm29. 
Thus, for the weak mixing, 9 is written as 

{A\H\B) _ {A\H,.,\B) 
^- AM ~ AM ■ ^^-^^ 

where iJint is a part of H which mixes the states \A) and \B). When we consider 

baryon mixings, 9 acquires spinor structures as will be discussed in the next subsec- 
tion. 

8 



3.2 Phenomenological ansatz 

Let us consider the case that A and B are baryons with different ffavor quantum- 
numbers and different masses. We further assume that the interaction which induces 
the mixing is small and can be treated in the 1st order perturbation. Then the 
correlation function Eq. (|2.1|) near the mass-shell of baryons A and B may be written 
as 

n^ (p I g) = AaAb -(eAM) -. (3.6) 

where we omit the single pole contributions. Here \a is a coupling strength defined 
by (0|?7a(0)|v4(p, s)) = AaMa(p, s), where ma(p, r) ma(p, s) = 2MA5rs and \A{p,s)) 
satisfies the covariant normalization. The phase of the state \A{p^s)) (|-B(p, s))) is 
chosen so that A^ (Ab) becomes real. Aa, Ab, Ma, Mb and AM = Mb — Ma take 
their unperturbed value in the 1st order of the mixing parameter G. A diagrammatic 
illustration of Eq. (|3.6| ) is shown in Fig.|l|. The state vector |g), which will be later 
identified with the isospin-asymmetric nuclear medium, is the major physical source 
of the mixing. In general, B is a 4 x 4 matrix in the spinor space. According to the 
discussion in Sec.^, the correlation function ( |3.6| ) must have a form 



n^ (p I g) = n| + n^V + n^'g + U^'{ia^,p^qn, (3.7) 

where Hip = Hip {p'^,p-q, q^) and 



n^ 



A^B 



n^ (/ = 5, \/i, V2) , n^\ = -n^^ (3.8) 



Ti 



A^B 



Thus the Eq.( |3.7| ) restricts the structure of 0. In fact, it cannot contain 7^, 7^7^ and 



0"^^, and one obtains 



= ^^ + e^'p + e^^q + e^'iia^.p^q"), (3.9) 



where the parameters 6^ {I = S,Vi,V2,Ti) are real. 

Now, let us consider an effective Lagrangian £int which describes Eq.(|3] 



A„t(x) = (^a(x) e ^b{x) + ^b(x) ^a(x)) am. (3.10) 

Here = 7°0^7° and \1'a (^s) is the field that describes the particle A{B). By 
the parity and time-reversal invariance of £int, ^a and \E'b have the same trans- 
formation matrices P and T under the parity {V) and time-reversal (T) transfor- 
mations. This is the same constraint for the interpolating operators tja and tjb 
in Sec. p.l| . Therefore, the fields \E'a and \1/b can simultaneously satisfy the rela- 
tions (0|^a(0)| A(p, s)) = ua{p,s) and (0 |^ij(0)| S(p, s)) = ub{p,s) with the states 



9 



\a G am \b 
Va X * * X riB 



A 



B 



State with q^ 



Fig.l: A diagrammatic illustration of Eq. (|3.6| ). The baryon A{B) couples to the 
interpolating operator tja (jIb) with the coupling strength Xa (As). The baryon mixing 
AM is induced by the state with four-momentum g'^ which will be later identified 
with the isospin-asymmetric nuclear medium. 

\A{p, s)) and \B{p, s)) defined below Eq.(|;6D. Then, {B \Hi^t\ A) = {A \Hi^t\ B) and 
(B\H,^,\A) = (^\H,^,\B)\io\d. 

We have implicitly assumed that the interaction does not contain any deriva- 
tives, which is equivalent to the assumption that O depends only on g^ and not on 
p^. Then, reduces to a simple form, 

Q = e^ + e^-'q (3.11) 

For later convenience, we define a dimensionless parameter 6^ as 

e^ = e^''q\ (3.12) 

We can make physical interpretation of the mixing angles Q^ and Q^ as follows. 
Consider the particle at rest in the rest frame of the medium (p = q = 0). Then the 
off-diagonal matrix elements of the Hamiltonian H\^x. (= — \ (^^ l^i-ai^^'x)) read 



{A\E,^,\B) = -[9'' + 9'') AM, (3.13) 

A \Hint\ b) = - (9^ - 9^) AM. (3.14) 



where 1^4) and \B) are normalized as (A|A) = {B\B) = 1. Together with Eq.( |3.5| ), 
we thus find that 9 (9) defined below corresponds to the mixing angle in the particle 
(anti-particle) channel: 

9 = 9^ + 9^, (3.15) 

9 = 9^-9^. (3.16) 

9^ and 9^ have formal analogy with the scalar and vector self-energies (S'^ and S^) of 
the nucleon in the nuclear medium. The nucleon and the anti-nucleon feel an optical 
potential S"^ + S^ and S"^ — E^ respectively p!8 |. 



10 



Let us rewrite the correlation functions in terms of the mixing angles 6^ and 
6^ defined above. Since Rell^ has a form deduced from Eg . (13 .61) and Eq.( p.ll 
obtain 



we 



Renf 



Ren^i 



Ren^^ 



ReUJ> 



XaX 



A'^B 



e^MA + e^^{j>-q) 9^MB + 9^^{p-q) 



XaXb 0' + r^ 



{p-qY 



Ml-p'^ 
1 1 



\a\b ( e""'' 



AM\ 



A^Ab -r^ 



2M) 
AM\ 



M J [Ml-p"^ Ml-p"^ 
Ma ^ Mb 



M|-p2 Ml-p"^ 
1 1 



2M 



} 



Ml-p2 Ml-p"^ 



(3.17) 
(3.18) 
(3.19) 
(3.20) 



We mention here that Relly given above satisfy Eq.( p.8D , since the mixing angles 
have the property 



e 



sy 



A^B 



-^^'^, 



which is confirmed by Eqs. (|3.13|) and (|3.14|) . 

For p^ = (cj, 0) and q^ = (g°, 0), the correlation function ( |3.7D is simplified to 



V I. ,\ „o 



n^. (u;) = n^ (ct;) + n^ (tu) i 



and its real parts become 



[tM = 


- A^Ab 




9^Mb + e^u 


^M = 


= A^Ab 




O^u + O^Mb 
M|-w2 



(3.21) 

(3.22) 
(3.23) 



UJ 



We decompose Uj^ {uj) further into even and odd parts 

I = S,V. Then we obtain a simple formula to be used later 

Ma M, 



Ren^^""^ (s) 



Ren^^°^ (s) 



Xa^b 
XaXb d 



B 



Ml- s Ml- s 



M\- s Ml- s 



(3.24) 

(3.25) 
(3.26) 



where s = uj"^ and I = S,V . 
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4 OPERATOR PRODUCT EXPANSION 

In this section, we carry out the operator product expansion (OPE) of IIj. up to 
dimension 4 and evaluate the in-medium matrix elements of local operators. As is 
well-known, the operators with Lorentz indices should be retained since they do not 



vanish in the medium [|T4|. Furthermore, we need to keep not only the iso-scalar 
operators but also the iso-vector ones to take into account the isospin asymmetry in 
the medium. 

4.1 OPE for S° — A mixed correlation function 

Taking A = A and i? = S", the retarded correlation function ( ^.2| ) reads 

U^{p\q) = t Jd'xe'P^q\R[r]j,{x)r]^om\q) , (4.1) 

which satisfies the dispersion relations (p.33|) , (|2.35|) and ( |2.36|) . On the other hand. 



the time-ordered correlation function ( |2.1| ) reads 

U^{p\q) = z Jd'xe'^^qn^{x)r]^om\q) , (4.2) 

which is useful for making OPE. \q) is the state vector corresponding to the nuclear 
medium with total four-momentum q^. 

For the interpolating operators r]\{x) and riY;o(x), we adopt the loffe's current 
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where ip{x) is the quark field with flavor ip, C denotes the charge conjugation matrix 
and a, 6, c are color indices. This current is symmetric under the exchange of ipi and 
ip2, i.e., T]^^i)2ip-i{^) = '7V'2 5/'i^3(^)- Thus, rj^{x) and rjY,o{x) may be written as 

riK{x)=isJ-{r]usd{x)-ridsu{x)} , rjj:o {x) = i V2 rj^dsix) ■ (4.4) 

These are the same interpolating operators used in the analysis of S" — A mixing in 
the vacuum [|l^, |l^. Under the time- reversal, ?7a(2;) and ?7so(a;) transform in the same 
way as the quark field ipi^)- 

Using Eq.( |4.4|) and the above mentioned exchange property, the mixed correla- 
tion (14.21) becomes 



2 

Ht (p k) = -yl 



d'xe'^' {q\T[r^usd{x)Vudsm\(l) - {n ^ d} 



(4.5) 



We carry out OPE of this correlation up to dimension 4 operators. The quark masses 
are kept up to the 1st order. Since Eq.( [4.5| ) is anti-symmetric under the exchange of 
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u and d, u — d symmetric terms such as the gluon condensate ( — G^ \ do not appear. 
All the diagrams contributing up to the order we consider are drawn in Fig.0. The 
explicit forms of OPE may be summarized as follows: 

n^ {p\q)=U^ + Yil'p + Yil'q + III' {la^^p^q'). (4.6) 



Ren| = ^ 



U-n'' 



{rrid - niu) p'^ log{-p^ 



1 /^ 



iTT^ 



dd — uu) p log(— p 



- —^ \ (2m„+mrf) {uviu) 



Ren^i 



2 

71 



1 



247r2 
1 



+ 



1 

ri-n' 



mu + 2md) (dvidj + {md-niu) {si/is) \ {n-p) log(-p^) 
dd — uu) {n-p) log(— p^) 
2 , nis (dd - uu^ + (m^ - m„) (ss) [ log(-p^) 
- <^ {nid (dd\ - rriu {uu) \ 



— 4 /d{n-iD)fid — u{n-iD)fiu) > log(— p^) + 



2{n-pf 

p2 






487r2 



dfid — uyiu) p^ log(— p^) 



— -^ I {nid (dd^ - rriu {uu) 
— 4 (d{n-iD)ftd — u{n-iD)ftu 



{n-p) log(— p^ 



Ren^ 



2 



iic^ 



{nig — nid) {uyiu) — (m^ — niu) (dyid 



+ {nid - rriu) {syis) \ log(-p^) 



.(4.7) 



(4.8) 



Nq' 



Nq\ (4.9) 



(4.10) 



where we have defined a normal vector n^ = q^/\/q^ characterizing the nuclear 
medium. Also we have replaced the in-medium matrix elements {q\- ■ -{q) by (• • ■) for 
simplicity. 








Fig .2: Diagrams of OPE up to dimension 4 for the S" — A mixied correlation function, 
(m stands for the u, d, s quark masses and V fo^' the covariant derivative.) 

The in-medium expectation values of isospin-asymmetric operators beyond di- 
mension 4 have large uncertainties. For example, (d\a-G)d — u\a-G)u\ which may 
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give a major contribution to the OPE at dimension 5, is not known. In fact, even its 
isospin-symmetric partner (q'^{cr-G)q) has large error (— 0.33GeV^ ~ +0.66GeV^)-/>N 
with p,^ being total nuclear medium density |T^. Therefore, we limit ourselves to the 
isospin-asymmetric operators up to dimension 4 in this paper. 

In the rest frame of the medium with n^ = (1,0), the decompositions to even 



s,v 



and odd parts of Refl-p' are written as 



Ren 



S(E) 



Ren, 



y(E) 



s, IpI) 



\s,\p\) 



1 





2 


IGtt 


v^ 


1 




2 


167r 


v^ 



1 (-471"^) 

-5mp'^log{-p^) H {5iqq))p^\og{-p'^ 



(4 



n 



Ren^^°) (s, IpI) 



Ren, 



V{0) 



A IpI) 



levrV ^/3 



2 

71 



4 
(4vr^) 



-7^ —^{^{qq))p iog(-p ) 



, (4.11) 
(4.12) 



5m Uq'^q) — (s^s)) — rn (S{q^q)) ( log(— p 



IGtt^ 



(4 



TT 



[ms {S{qq)) + 5m (ss) + X) log(— p^ 



-Ha- I 



pz 



(4.13) 



(4.14) 



with 



X ^ -{{6{qzPq)) - A{6{qzDYq))} 



and p^ = s — p^. In the above formulas. 



(4.15) 
denotes the u — d average, namely 



m = {m,u + mid) /2, q'^q = \y)u + Sd\ /2 and qq = [uu + dd\ /2, while 5 denotes 
d — u difference, namely 5m = m^ — m^j, 5{qq) = dd — uu and (5(g^g) = (i^d — u'''u. 

4.2 In-medium condensates 

In the previous subsection, we have encountered various u ~ d symmetric and u — d 
anti-symmetric condensates. In this subsection, we will evaluate those in a model 
independent way using the low density expansion. 

First of all, the expectation value of the local operator O has a vacuum part 
which is density independent and the medium part which is density dependent; 

{0) = {0)^ + {0)^^. (4.16) 

At low density, {O)^^^ is expanded as 

{0)^. = {0iPp + {0lPn + ---, (4.17) 

where Pp (pn) is the proton (neutron) density, and {O) {{O),^) is the spin-averaged 
expectation value taken by the one particle state of the proton (neutron). 

{0)j, = Jdhm \0\ N) - (0 \0\ 0)} , (4.18) 
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where (A^| A^) = 1 for N = p,n. 

Since the vector condensate (q^q) is nothing but the quark number density, we 
have 

(?g) = ^Piv, {6iqk)) = Sp^, {sh) = 0, (4.19) 

where pr^ (= Pp+Pn) is the total nucleon density, 6pj^ (= Pn—Pp) is the n—p asymmetry. 
The scalar condensate {qq) = ijlq)^ + {qq)Mcd i^ evaluated by the Feynman- 
Hellmann theorem |]23[| ; 

\ / Med 

where A is a parameter in the QCD Hamiltonian, and S is the energy density of the 
nuclear medium. At low density, 

S = (t^QCD) c^MpPp + Mn Pn = MnPn + SMn 6p^/2 (4.21) 

\ / Med 

where M^ = (Mp + M„)/2 and 6M]m = M„ - Mp. Under the choice X = 2m together 
with the mass term of the QCD Hamiltonian 

^Sas? = rn{uu + dd) + —{dd - uu) + tUs'ss, (4.22) 



one finds from Eq.( |4.20| ) 






where we have used a definition of the nucleon sigma-term a^ = m 



dm 



On the other hand, the choice A = 5m/2, with Eq. ([4.20|) gives 

which is valid up to the 1st order in 5m. 

The strange-quark condensate up to the 1st order in p^v reads 

(^s)Med - \ ((^s)p + (^^)n) Pn = \ [{vq\ + {fq\) Pn. (4.25) 

Here y is a parameter characterizing the OZI violation in the nucleon 0. 

Thus the quark condensates (qq) in isospin-asymmetric nuclear medium are 
summarized as follows: 

mq)) - <^(^^))o + ^^^P- (4-27) 

{SS) ^ {ss)^ + y^p,. (4.28) 
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By using the result of evaluation of (S^qiD'^Y^q)) in Appendix, we obtain 



X ^ -{{6{qtPq))-4{6{qzDYq))} 



(4.29) 



where AKp"^) is the 2nd moment of the parton distribution function of the proton 
defined by Eq.( |A.4|) in Appendix. 



4.3 Summary of OPE for Rellrp 

In the isospin-asymmetric nuclear medium, the OPE of the correlation function is 
finally expressed as 



Ren: 



5(E) 



s, IpI) 



/ 1 



Ren^^'') (s, \p\) 



VievrV V3 
/ 1 \ 2 



1 (-4vr^) 



VievrV Vs 



(4 



n 



6{q^q))p^\og{-p^ 



Ren^p (s, IpI) = 0{m6pis,, 6mpi^), 



(4.30) 

(4.31) 
(4.32) 



Ren, 



V(0) 



's, \p\) 



(47r2 



>s {^iQQ))o + ^^ (s^)o) log(-P^ 



167rV V3 

+ ^ S{Al{p^)) (Mj, 6p^ + SMn p,/2) llogi 
+ 0{ms5pt,, 6m pr,), 



-P 



2 A 

3 p2 



(4.33) 



where p"^ = s — p^. 

Above formulas are valid up to the 1st order in the quark masses and the baryon 
density. In other words, the terms such as 0{6m) and 0{5pm) are kept, while the 
terms such as 0{mu^d,sPN), 0{mu,d,s6pN) are neglected. Within this approximation, 
the scalar-odd correlation n^- vanishes as shown above and cannot be used to 
construct sum rules. Since it is theoretically consistent to use the correlation functions 
with the same reflection symmetry under u <-> —u, we will exclusively use the "even" 
correlations n^. and n^- in the following analyses. 



5 QCD SUM RULES 

In this section, we construct finite energy sum rules (FESR) |T^ and Borel sum rules 
(BSR) [|] on the basis of the retarded correlation function (|4.1| ) and the dispersion 
relations (|2.35|) and (|2.36| ). For the phenomenological side, we use the ansatz given 
in Eqs.QOSl) and (^^. The OPE side is given in Eqs.(|Al]-|l§. 
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5.1 Finite Energy Sum Rules (FESR) 

In FESR, we identify the integral of Imll^ extracted from OPE with that introduced 
phenomenologically: 

„i(0) „i(0) 

J^ " ds {imnS {s, \p\) V~s} ^" = y/ ds {imnS (s, \p\) V^} ^", (5.1) 

UE) 1(E) 

r ds{lmU'iSnis,\p\)/V~s}s- = r ds{lmU'gl{s,\p\)/V~s}s-, (5.2) 
Jo '^ -' Jo '- -' 

where s = uj'^. So and 5*0 are the continuum thresholds of Imllpjjgjj and Imllpj^gj^ 
respectively. These phenomenological spectral functions at p = in the left hand 
side of the sum rule are obtained as follows. We substitute the ansatz ( |3.25| ) and 
( lOGi) with A = A, B = J:° into Eqs.(^) and (lOSi ) to obtain 



Imll 



Phcn 



imn 



V 
Phen 



[UJ 



[UJ 



— AaAso 



TC 



-AaAso 



e"" +e^j{6{uj-MA)-6{uj- Ms..)} 
- (e^ -e^){6{uj + Ma) -6{uj + Mso)} 
[e^ + e^){6{uj- Ma) -6{uj~ Meo)} 
+ (^^ -9^){6{u + Ma) -6iu + Mj,o)} 



(5.3) 



(5.4) 
The even-odd decompositions of the above formula give Imllpj^gj^ and Imflpj^j^. As 



we have discussed in Sec.2]2| and Sec.[4.3|, we consider only the first-order effect of 
the isospin-asymmetry on the mixed correlation function and neglect the effects of 
0{mu4,s Pn) and Oirriu^^s^PN)- Therefore, the pole positions Ma and M^o take their 
vacuum value and only the pole residues are affected in a different way in even and 
odd spectral functions. The situation is the also same for excited states of A and S°. 
Therefore, in the present approximation, the continuum threshold also takes their 
vacuum value (5*0 = 5*0 = Sq ) and only the height of the continuum is affected 
by the isospin asymmetry. 

The OPE motivated spectral functions at p = in the right hand side of 
the sum rule are obtained as follows. We substitute Eqs.( |4.1ll - [4.14| ) together with 
Eq. dO^) into Eqs.fglSTl) and (lO^ ) to obtain, 

1 \ 2 



Imn 



OPE 



\UJ 



— vr sgnu^ 



levrV V3 



mq))^' 



+ 



(47r2 



\5ni Uq^q 



— is' s) ] — m 






■ u 



ImHopE (uj) 



-n sgn(u;j 



levr^ 



2 
73 






(5.5) 



+ 



^ I (m, {6{qq)) + 6m {ss)) + ]- 5{Al{^Ji'')) [m^ Sp^ + 6Mm 9^/2) \ uj 



(5.6) 
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Here, sgn(ci;) = oj /\uj\ {u ^ 0) and sgn(O) = 0. The even-odd decompositions of the 
above formula give Imllop^ and ImllQpp,. The isospin-asymmetric condensates affect 
the magnitude (height) of even and odd components in a different way. 

As we have mentioned at the end of Sec.^, we use Rell^^ for the actual analysis, 
which corresponds to adopt the sum rule ( |5.1|) . Remember that the even part Rellj^ 
is related to the odd part Imllj^ through the dispersion relation (|2.35|) . Resulting 
FESR for p = reads 






2 
73 

2 
73 



1 Cn+S 

-5m ^^^ 
4 n + 3 



-47r2 



(4 



TT 



6{qh) 



Qn+2 

n + 2 



mq)) 



Qn+2 

n + 2 



(5.7) 
(5.8) 



where /3aso = IGtt'^AaAso. As is evident from the right hand side of Eq.( 



e' 



appears only in the nuclear medium, 6 — ^^^^. 



vacuum part and in-medium part 6^ 
from the sum rule (15.71) to obtain 



el 

35 



^0 + ^ivlcd- 



On the other hand, 9^ has both 
We subtract out the vacuum part 



/3aeo (mIU+i 



^j2n+l 



9 



2 
73 



-47r2 



mQ)i 



n + 2 



(5.9) 



Combining Eqs.(|5.8|) and (|5.9|), one has a simple formula for the ratio 9l^/9- 



s . 

Med" 



6{qk) 



9^ 



2 mq)). 



-2 



d{6Mj, 
d{5m) 



(5.10) 



Mn-Mp 



where we have used Eqs.( [4.11J| ) and ( |4.24D for the last equality with ^Mat 
and 5m = m^ — m^j. 

Let us estimate the right hand side of Eq. (|5.10| ). The n — p mass difference 
5Mn = 1.29 MeV is known to be decomposed into two parts, 5Mn ^ 6M§^ + 6M^^^ . 
Here (5M™(~ —0.76 MeV) originates from the electromagnetic interaction of 0{a), 
while SM^ (~ 2.04 MeV) is due to the u — d quark mass difference of 0{Sm.) 0. 
Therefore, in the leading order of Sm and a, one finds 



d{6M, 



N 



d{5m) 



d{5m) 



6m 



0,52 1^^) 

om, J 



(5.11) 



where we have used 5m, = 3.9 MeV as a typical value at the renormalization /i^ = 
1 GeV^ 1^, |12| (See also Table |I]). Thus we find that the scalar and vector mixing 
angles induced by the nuclear medium have opposite sign and approximately equal 
in magnitude. 



^Mod/^i; 



-1. 



(5.12) 



18 



This together with the definition of the total mixing angles Eqs.( p.l5| ) and ( |3.16| ) 
implies that the medium modification of the particle mixing is largely cancelled be- 
tween the scalar and vector, while the anti-particle mixing is enhanced in medium. 
The magnitude of the mixing angles will be discussed in SecB. 



5.2 Borel sum rules (BSR) 



In BSR, we make a Borel transform of the dispersion relations ( p.35|) and ( p.36|) for 



the retarded correlation (|4.1j ) in the deep Euclidian region s = u^ 



— oo: 



B 



IpOs'^nKUsMpDV?' 



B 



n Jo s' — s 

1 p^-^^, Imngg,(3Mp|)/v/? ' 
71 Jo s' — s 



where the Borel transform B is defined as 



-s 



B 



B 



d 



RenS(s,|p|) 
RenS(^,|p|) 



(5.13) 
(5.14) 



^P(^)]-4^oo(n-l)!V^. 



n(.). 



(5.15) 



with M being the Borel mass. 

The left hand side of the sum rule (the phenomenological side) is assumed to 
have the pole + continuum structure: 



A (E,0) 
■■Phen 



iipi'r^ (s, ipI) 



n' 



Z(E,0) 



/(E,0) 



'■Phcn(polo) \^i \P\) + nphek(cont) l'^' \P\) i 

where the continuum part Ilpj^g^/^^^^-j is extracted from ImllQp^ , 

im4SLt) (^, Ipl) = imnS^^ (., Ipl) e{.s - Sl^^'^^ 



(5.16) 



(5.17) 

with 6{x) being the step function. 5*0 and 5*0 are the continuum thresholds of 
Imllpjjjj^ and Imllpj^gJ^ respectively. The pole part for p = has been already discussed 
in Eqs.( |5.3| ) and ( p.4| ). The right hand side of the sum rule (the OPE side) is derived 
from Eqs.( ^4.1l| - [4.14|) . For the reason which we denoted in the previous subsection, 
we use universal threshold S^ = So —So . However the height of the continuum 
is affected by the isospin-asymmetric nuclear medium, and Imllpj^gj^. ^^^n (cu) becomes 
asymmetric under u ^->- —uj. This is taken into account in Eq.( |5.17D through the 
asymmetry of ImllopE {uj). Its explicit form is shown in Eqs. (|5.5| - |5.6|) . 

For the scalar mixing angle, resulting BSRs from Eq.( p.l3D at p = become 



/3aso (Mso 



g-A^^o/A^' 



Mac 



-Ml/M'- 



e' = f[M\So 



(5.1^ 
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with 

and for vector mixing angle 

/3aso (Mso 

with 



2 1 VmV 



f^\ ^ (- 



47r^)(5(gg))M^Ei(A) 



M^JM^ 



/"K.«")-^ 



M^e-'-'i'"')0" ^f{M\S, 



(ip (,„.,)) M'^,(^) 



(5.19) 



Here E„(x) = 1 - e-^Er=oa;7^!- 

One can simphfy the above sum rules without loss of generality by expand- 
ing the left hand side in terms of a small parameter AM/M = (M^o — M/C)/M ~ 
(Mso - Ma)/((Mso + Ma)/2) ~ 77.0MeV/1.15GeV ~ 0.07. Using this expansion 
and neglecting 0{{AMY) contribution, Eq.( p.l8 ) becomes 



(m^ - 2M') ef = f^[m'^, So 



(5.20) 



where F^{M^,So) = f^{M^, So) M^e^'/^^' /{Pa^o AM) and we put the suffix 1 to 
the mixing angle for later convenience. One may alternatively take derivative of 
Eq. (|5.18|) with respect to M^ to enhance the lower dimensional operator in OPE and 
then expand the result by AM/M. Then we obtain the second sum rule for the 
mixing angle 

(SM^ - 2M') ^fi = F^[m\ So) , (5.21) 

where F^(M2, So) = {df{M\ So)/dM^){MyM' ) e*^'/^^V(/5Aso AM). 

Extracting a term which is proportional to M^ and a term which is M^- 
independent in Eqs.( p.20| ) and (|5.21|) , one finally arrive at four types of sum rules, 
which we call Type Ia,lb,IIa and lib. 

d 



nS _ 

^" rfM2 



f'{m\So) 



9i = {M'-^^F%M',So)-F'{M\So)]/{2M^ 



and 



^-4^^i^'^«) 



eL = \m'^^F'{m^So) - F^M^So)] /i2M^ 



(5.22) 
(5.23) 
(5.24) 

(5.25) 



The Borel sum rules for the mixing angle 6^ can be formulated exactly in the 
same manner starting from Eq. (|5.19|) . The results are obtained by replacing the suffix 
ShyV in Eqs.(p:22D, ( [OBI) , (jO^ ) and ( ^1251 ). 
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6 NUMERICAL RESULTS 

In this section, we evaluate the absolute value of the mixing angle with the use of 
the BSRs, Type la (lO^) , lb (g^, Ila (^M) and lib {^^, supplemented with the 



FESRs ( |5.7| ) and ( [5.8| ). To extract the mixing angles from the sum rules, we need to 
know various QCD parameters (vacuum condensates and quark masses), the coupling 
strength |/3aso| and also the Borel window in which Borel analysis is made. They are 
determined by the following procedures. 

1. QCD parameters and /^a^o in the vacuum 

In Table |l|, the QCD parameters which we use in our analysis are summarized. 
These parameters reproduce the mass spectrum of octet baryons in QCD sum 
rules within 10% [T^. (OPE up to dimension 7 and quark masses up to the 



2nd order have been taken into account in this analysis.) /3aso(= IQit^X^X-^o) 
has been determined by the BSR in the vacuum for "diagonal" correlations 
[A = B = A and A = B = S°). Using the parameters in Table |I|, we obtain 
|/3aso| = 2.5 GeV^ from the scalar-even sum rule [l^. (Note that the sum rules 



for the diagonal correlations provide only the absolute values of Aa and A^o.) 
The optimum threshold 5*0 turns out to be 3.2 GeV^ from the Borel stability. 
This number is consistent with the position of the second resonances of A and 
E°. 

2. Mixing angle in the vacuum 

The scalar-even BSR for S"— A mixing angle in the vacuum 6^ (in which OPE up 
to dimension 7 and quark mass up to 2nd order have been taken into account) 
provides l^ol = 1-4 x 10^^ (^r^) with the Borel window 1.4 GeV^ < M^ < 
2.6 GeV^ and the threshold So ~ 3.2 GeV^ ^. Zhu et al. ^ obtained |0o| ^ 
7 X 10^'^ for different QCD parameters with |/5aso| = 1-76 GeV^. We will discuss 
the effect of this difference to the in-medium mixing angle at the end of this 
section. 

3. QCD parameters and Borel window for in- medium mixing angle 

The QCD parameter essential for obtaining the in-medium mixing angles is 
Eq.O); 

d{SMN)/d{Sm) = 0.52. 

In the standard Borel analysis, the Borel window is chosen such that the higher 
orders in OPE and the continuum contribution are well suppressed. The Borel 
window satisfying these conditions for the S° — A mixing angle in the vacuum 
is 1.4 GeV^ < M^ < 2.6 GeV^ as shown above where OPE up to dimension 7 
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has been taken into account [|T^]. For the mixing angle in the medium, we have 
OPE only up to dimension 4 in which the medium effects appear only in the 
highest dimensional operators and are dominant. Therefore, it is difficult to find 
the Borel window and to reach similar level of Borel stability. This is shown 
in Fig.|^ where in-medium mixing angle Ofj^ as a function of the Borel mass for 
different values of 5*0 is plotted. Since the Borel curve is not enough stable in 
the medium, we simply adopt the Borel window determined in the vacuum and 
extract the in-medium mixing angle by making average over the Borel window. 



m (= (m„ + md)/2) 


7.7 MeV 


5m/ m {5m = m^ — mu) 


0.51 


ms/m 


19 


{ml{={uu + -dd)j2) 


(-275MeV)3 


^{={dd)j{uul-l) 


-5.6 X 10-3 


n={-ss)j{uu\-i) 


-0.22 


{^G\ 


(347 MeV)^ 


ml {= -9s{qi(r-G)q)J {qql) 


0.91 GeV^ 



Table 1: QCD parameters in the vacuum at the renormalization scale 1 GeV^. Those 



are determined to reproduce the octet baryon spectrum [12 
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Fig. 3: The scalar angle in the Type Ila sum rule is shown as a function of the Borel 
mass M^ for different values of the threshold So [GeV^] . The straight horizontal lines 
imply the Borel window and the averaged value in the Window. Pah^ = 2.5 GeV^ 
and 5m, = 3. 9 MeV are used. 
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Adopting the procedure described above, we obtain the in-medium mixing an- 
gles 9^'^ as a function of the continuum threshold 5*0 for four different types of BSR; 
Type la, lb, Ila and lib. They are shown in Fig.^ and in Fig.^ at nuclear saturation 
density p^ = 0.16 fm~ = po and at typical value of the isospin-asymmetry for heavy 
nuclei such as Pb, Spff/p^f = 0.21. 
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Fig.4: The scalar angle is shown as a function of the continuum threshold Sq [GeV^] . 
The left panel is for BSR Type I and the right panel is for BSR Type II with corre- 
sponding n-th order FESR. /3aso = 2.5 GeV^ and 5m = 3.9MeV are used. 
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Fig.5: The vector angle is shown as a function of the continuum threshold 5*0 [GeV^]. 
The left panel is for BSR Type I and the right panel is for BSR Type II with corre- 
sponding n-th order FESR. /?aeo = 2.5 GeV^ and Sm = 3.9MeV are used. 

To reduce the uncertainties due to the absence of higher dimensional operators 
in OPE, we examine the reliability of each type of sum rules in the following ways. 
First of all, if the BSRs are consistent with each other, 9^'^ as a function of 5*0 should 
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have a similar behavior between Type la and Type lb and between Type Ila and 
Type lib. Such comparison is shown in Fig.|^ and Fig.^. The figures indicate that 
Type II sum rules are more reliable than Type I from this criterion. 

One can make further selection of a reliable BSR by the comparison with FESR. 
Remember that the n-th order term in the 1/M^-expansion of the BSR is equivalent 
to the n-th order FESR. Therefore, if OPE is well behaved, BSR and corresponding 
FESR should give the same result. Such comparison is also shown in Fig.^ and 
Fig.|^. From the right panels of Fig.| and Fig.^, we conclude that the Type Ila is 
more reliable than Type lib for reasonable range of the continuum threshold located 
around the second resonances of A and S°, 5'o — 3.2 GeV^. 

In Fig.^ we show the scalar angle 6^ and the vector angle 6^ in the Type 
Ila as a function of the continuum threthold 5*0. The curves in the Fig.|^ indicate 
the maximum (Max), the minimum (Min) and the average (Avg) value in the Borel 
window. Fig.|^ shows the n — p asymmetry dependence of 6^ and 6^ for the total 
density pr^ = 0.5 po, po and 1.5 po. 
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Fig.6: The scalar and vector angles are shown as a function of the continuum thresh- 
old 5*0 [GeV^] for Type Ila. The curve with the label "Avg" is the average of the sum 
rule over the Borel window. The curve "Max" ( "Min" ) is the maximum (minimum) 
value of the sum rule in the Borel window. 



Finally, by using the Type Ila sum rule, we obtain the scalar and vector mixing 
angles as follows; 



9 



[-(0.19 ± 0.02) a„p(p^/po)] 
[+(0.20 ± 0.02) a„p(p^/po)] 



2.5 GeV^ 

'2.5GeV^^ 

, /5aeo 



'3.9MeV\ 
5m J 



(6.1) 
(6.2) 



24 



where the n—p asymmetry is defined as a„p = Sp^/ Pn- The error bars are determined 
by the minimum and maximum values of the sum rule for mixing angles in the Borel 
window. For typical values of the parameters, 5m = 3.9 MeV and |/3aso| = 2.5 GeV^, 
Pn = Po and anp = 0.21, 6*^^^ dominates over the vacuum mixing angle Oq. Also, the 
relation OY^^^/O^^^^ ~ —1 discussed in Sec.| is well satisfied. 
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(^-nv (= 5pi^/ Pf^) for different values of the total density pjv- 



The particle and anti-particle mixing angles 9 and 6 are obtained from Eq. ( p.l5| ) 
and ( p.l6|) as 



9 = [(0.01 ± 0.04) a„p(p^/po)]sgn(/3Aso) + ^o, 
9 = [(-0.39 ± 0.04) a„p(p^/po)]sgn(/3Aso) + eo, 



(6.3) 
(6.4) 



with 5m = 3.9 MeV, |/3aeo| = 2.5 GeV^ 

In the analysis of the S° — A mixing in the vacuum by Zhu et al. [Ill], they 
use different set of QCD parameters from Table ^ (in particular 5m, = 3.0 MeV 



and (qq 



-241 MeVf) and obtain |/?, 



'AEO| 



1.76 GeV*^ with the threshold 5*0 



3.4 GeV^ and the Borel window 1.3 GeV^ < M^ < 2.5 GeV^ Substituting these values 
into Eqs.(|6.1D and (|6.2|), one obtains 



9 = [(-0.12 ± 0.08) «„p(p^/p„)]sgn(/?Aso) + ^o, 
9 = [(-0.75 ±0.09) a„p(p^/p„)] sgn(/3Aso) + ^o. 



(6.5) 
(6.6) 



which are qualitatively consistent with the result obtained using our parameter set. 
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7 CONCLUSIONS 

In this paper, we have studied the S" — A mixing angles in the isospin-asymmetric 
nuclear medium by using the QCD sum rules. 

Firstly, we have discussed general properties of diagonal and off-diagonal cor- 
relation functions of the baryonic currents. We found that the off-diagonal (mixed) 
correlation function consists of scalar, vector and tensor terms. They are further de- 
composed into even and odd parts in terms of the reflection symmetry under uj <-^ —uj. 
Then we derived dispersion relations for each component. 

Secondly, we examined the general structure of the mixing angle for baryons and 
introduced two independent mixing parameters 6^ and 6^ for the baryon at rest inside 
the medium. The sum (difference) of these parameters are shown to be the particle 
mixing angle 9 (the anti-particles mixing angle 9). This situation is analogous to the 
self-energy of the nucleon and anti-nucleon in the relativistic mean-field theories. 

Thirdly, we have carried out the OPE for the S" — A mixed correlation function. 
Then we constructed sum rules for 9^ and 9^ . From the finite energy sum rules, 
we found that ^'J^^j/^Mcd ~ ~^- This implies that the particle mixing angle 9 {= 
qS _|_ qV^ jj^ ^]^g medium is nearly equal to the one in the vacuum, and the isospin- 
asymmetric medium affects mainly the anti-particle mixing 9 {= 9^ — 9^). From the 
Borel sum rules, we evaluated the in-medium parts of 9 and 9 numerically. The results 
are summarized in Eqs. (|6.3| - |6.6|) in Sec.^ As the baryon density and the isospin- 
asymmetry of the medium increase, the anti-particle mixing is enhanced, while the 
particle mixing remains less than 20% of the anti-particle mixing. 

The strong correlation between 9^^^ and 6'^^^ and the strong modification of 
the anti-particle mixing in the isospin-asymmetric medium shown in this paper are 
model independent consequence supported both by the finite energy sum rules and 
the Borel sum rules. On the other hand, the absolute magnitude of each mixing angle 
has uncertainties due to the absence of higher dimensional operators in OPE. Better 
evaluation of the matrix elements of isospin-asymmetric operators beyond dimension 
4 is necessary for precise determination of the mixing angles. Also, it is an open but 
interesting problem to study whether one can measure the anti-particle mixing in 
nuclei in the laboratory experiments. 
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APPENDIX 

In this Appendix, we evaluate a dim. 4 isospin anti-symmetric condensate {5{qiD°'^°q)) 
following |l^ . 



In the medium, the dim. 4 quark condensate (qiD^Y^q) is represented as 

(qiDYq) = \m, {qq) + {qt{DY- ^/4) q) ■ (A.l) 

The vacuum part of the second term vanishes 

(g^(DV-^/4)g) ^ (g^(DV-^/4)g)pPp+ (gzpV-^/4)g)„Pn. (A.2) 



The expectation value taken by the proton at rest is [jl4| 

mDY- ^/4) q), = I M.AUfx') (A.3) 

where AK/i^) is the 2nd moment of the parton distribution function q {x, fi'^), q {x, fi"^) 
of the proton 



Al{^^') = 2j dxx"-i{g(a;,/i2) + (-l)"-ig(a;,/i2)}. (A.4) 



rl 
/O 

For isospin anti-symmetric condensate {6{qiD°'j''q)), we obtain 

{5{qi{DY- ^/4) q)), = ^ M,5(Al{^^')). (A.5) 

Then, the expectation value taken by the neutron at rest becomes 

{5{qi{DY- ^/4) q)l = -^ MJ{Al{^^')) (A.6) 

up to the 1st order in isospin-asymmetry. Thus we finally arrive at 

{6{qzDYq)) ^ I {m{6{qq)) + Sm{fq)} -^6{Al{fi')) {M^ Sp^, + 6M^ p^/2). (A.7) 
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